The structure of Inter-Urban traffic: 
A weighted network analysis 
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We study the structure of the network representing the interurban commuting traffic of the 
Sardinia region, ftaly, which amounts to 375 municipalities and 1,600,000 inhabitants. We use 
a weighted network representation where vertices correspond to towns and the edges to the actual 
commuting flows among those. We characterize quantitatively both the topological and weighted 
properties of the resulting network. Interestingly, the statistical properties of commuting traffic 
exhibit complex features and non-trivial relations with the underlying topology. We character- 
ize quantitatively the traffic backbone among large cities and we give evidences for a very high 
heterogeneity of the commuter flows around large cities. We also discuss the interplay between the 
topological and dynamical properties of the network as well as their relation with socio-demographic 
variables such as population and monthly income. This analysis may be useful at various stages in 
environmental planning and provides analytical tools for a wide spectrum of applications ranging 
from impact evaluation to decision-making and planning support. 
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I. INTRODUCTION 



The informatics revolution has enabled the systematic 
gathering and handhng of a large array of large scale net- 
works data sets, enabling the analysis of their detailed 
structural features. In particular, mapping projects of 
the WWW and the physical Internet offered the first 
chance to study the topology and traffic of large-scale 
networks. Gradually other studies followed describing 
population networks of practical interest in social sci- 
ence, critical infrastructures and epidemiology (Albert 
and Barabasi, 2002; Newman, 2003; Dorogovtsev and 
Mendes, 2003; Pastor-Satorras and Vespignani, 2004). 
The possibility of accessing and mining large-scale data 
sets allows a more detailed statistical analysis and the- 
oretical characterization of correlation patterns, hierar- 
chies and community structure of these networks, along 
with the setup of new modeling frameworks. In particu- 
lar, the systematic study of these systems has shown the 
ubiquitous presence of complex features, mathematically 
encoded in statistical distributions with heavy-tails, di- 
verging fluctuations, self-organization and emerging phe- 
nomena and patterns. 

The statistical properties induced by these complex 
features indicate the presence of topological and struc- 
tural properties that do not find an explanation in the 
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paradigm put forward by Erdos and Renyi with the ran- 
dom graph model (Erdos and Renyi, 1960). Indeed, even 
if the Erdos-Renyi model rationalizes the small-world 
property of networks (i.e. the short distance measured 
in number of links among nodes), it fails in reproduc- 
ing the high level of local cohesiveness observed in many 
networks (Watts and Strogatz, 1998). Additionally, sev- 
eral of these networks are characterized by the statistical 
abundance of nodes with a very large degree fc; i.e. the 
number of connections to other nodes. For these 'scale- 
free' networks, the degree probability distribution P{k) 
spans over a wide range of k values which signals the ap- 
preciable occurrence of large degree nodes, the 'hubs' of 
the system. Finally, it is important to remark that the 
topological features of networks turn out to be extremely 
relevant since they have a strong impact in assessing their 
physical properties such as robustness, vulnerability, or 
the ability to spread a disease (Barabasi and Albert, 
1999; Dorogovtsev and Mendes, 2003; Pastor-Satorras 
and Vespignani, 2004). 

More recently, the activity on complex networks has 
been extended to the characterization of weighted net- 
works. This representation allows the consideration of 
features pertaining to the dynamics and traffic flows oc- 
curring on networks, adding another dimension in the 
description of these systems. Also in this case, the 
analysis and the characterization of weighted quanti- 
ties have pointed out the presence of large scale hetero- 
geneity and non-trivial correlations (Barrat et al, 2004). 
The weighted graph representation provides valuable el- 
ements that open the path to a series of questions of fun- 
damental importance in the understanding of networks. 
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Among those, the issue of how dynamics and structure 
affect each other and the impact of traffic flows on the 
basic properties of spreading and congestion phenomena. 
Finally, it spurs the more theoretical questions on how 
all these properties may be considered in generative net- 
works models. 

In the realm of urban and environmental planning, spa- 
tial analysis and regional science, the interest for com- 
plex networks has noticeably increased during the last 
years. In his investigation of the structure of the mod- 
ern urban fabric, Salingaros (2001 and 2003) invokes the 
concepts of small world network and scale free proper- 
ties and suggests a parallelism between the city and an 
ecosystem. Starting from the principle that connective 
webs are the main source of urban life, Salingaros invites 
planners to apply complex networks analysis to under- 
stand properties such as resilience and self-organization. 
While these studies are mostly based on conceptual ar- 
guments, other authors (Shiode and Batty, 2000; Batty, 
2001 and 2003) promote the development of quantitative 
analyses with a twofold perspective: (i) explaining the 
behavior of complex networks such as the world wide web 
in terms of geography, demography and economics; and 
(ii) improving the efficiency of advanced spatial analysis 
GIS-based tools, by integrating the foundation of com- 
plex networks analysis into geographic information sci- 
ence. Batty (2001) along with Salingaros acknowledges 
the existence of a link between efficient urban connectiv- 
ity and self-preserving behaviors and fosters the develop- 
ment of complex networks to study the functional mech- 
anisms of cities. It is interesting to note that a similar 
approach has been extremely fruitful in the study of the 
Internet in relation with the geographical and social en- 
vironment (Yook et al, 2002; Gorman, 2001; Barthelemy 
et al, 2003; Gorman and Kulkarni 2004; Pastor-Satorras 
and Vespignani, 2004; Schintler et al, forthcoming). 

In the domain of geographic information science, Jiang 
and Claramunt (2004) proposed to use complex network 
tools to investigate quantitatively the urban space syn- 
tax. Aiming at the improvement of the effectiveness 
of network GIS advanced analysis, they use topologi- 
cal measures of connectivity, average path length and 
generalised clustering to compare three large street net- 
works of the cities of Gavle, Munich and San Francisco. 
Fertile research directions are also currently explored in 
the domain of transportation systems and human traffic 
analysis. For example, the small-world properties of the 
Boston subway system have been characterized by Latora 
and Marchiori (2000-2003) and Sen et al (2003) studied 
the topology of the Indian railway network. At the urban 
level, a network study has been carried out by Chowell et 
al (2003) using large scale simulations (with census and 
demographic data integration) in order to describe urban 
movements of individuals in Portland, Oregon (USA). 
On a larger scale, recent studies have focused on the net- 
work properties of major transportation infrastructure 
such as the US Interstate highway network and the air- 
port network. (Gastner and Newman, 2004; Guimera et 



al, 2003). Finally, Barrat et al (2004) have provided the 
first study of the worldwide airports network including 
the traffic flows and their correlation with the topologi- 
cal structure. In this research area, the network approach 
might provide relevant information on issues such as traf- 
fic analysis and risk assessment in case of damages and 
attacks. 

In the present work we use a network approach to 
study the Sardinian inter-municipal commuting network 
(SMCN), which describes the everyday work and study- 
led movements among 375 municipalities in the Italian 
region of Sardinia. We obtain a weighted network rep- 
resentation in which the vertices correspond to the Sar- 
dinian municipalities and the valued edges to the amount 
of commuting traffic among them. We provide a detailed 
quantitative study of the resulting network, aiming at 
a first characterization of the structure of human traffic 
at the inter-city level and its relation with the topologi- 
cal structure defined by the connectivity pattern among 
cities. Very interestingly, while the network topology ap- 
pears to fit within the standard random graph paradigm, 
the weighted description offers a very different perspec- 
tive with complex statistical properties of the commuting 
fiows and highly non-trivial relations between weighted 
and topological properties. We discuss these properties 
in relation with the geographical, social and demographi- 
cal aspects of the system and show the potentiality of the 
network approach as a valuable tool at various stages of 
the policy-making and environmental planning process. 

The paper is organized as follows. In the next sec- 
tion, we describe the dataset. In section 3, we carry out 
the analysis of the topology of the network. In section 
4, we analyze the traffic and weight properties of this 
network and in section 5, we discuss the implications of 
our network analysis results in terms of urban and en- 
vironmental planning. In section 6, we compare topo- 
logical and weighted network analysis results with socio- 
economic characteristics of the Sardinian municipalities 
and we discuss their relations. We finally summarize our 
results in section 7, while in section 8 we discuss the dif- 
ferent research directions that this work suggests. 



II. SETTING THE CASE STUDY: DATA AND 
GEOGRAPHICAL FEATURES 

Sardinia is the second largest Mediterranean island 
with an area of approximately 24, 000 square kilometers 
and 1, 600, 000 inhabitants. Its geographical location and 
morphological characters have determined an important 
history of commercial and cultural relations with trans- 
borders external communities and have also favoured the 
development of a strong social identity and of important 
political movements toward self-reliance and autonomy. 
At the date of 1991, the island was partitioned in 375 
municipalities, the second simplest body in the Italian 
public administration, each one of those generally cor- 
responding to a major urban centre (in Figure 1, on the 
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left, we report the geographical distribution of the munic- 
ipalities). For the whole set of municipalities the Italian 
National Institute of Statistics (Istat, 1991a) has issued 
the origin-destination table (ODT) corresponding to the 
commuting traffic at the inter-city level. The ODT is 
constructed on the output of a survey about commut- 
ing behaviors of Sardinian citizens. This survey refers 
to the daily movement from the habitual residence (the 
origin) to the most frequent place for work or study (the 
destination): the data comprise both the transportation 
means used and the time usually spent for displacement. 
Hence, ODT data give access to the flows of people regu- 
larly commuting among the Sardinian municipalities. In 
particular we have considered the external flows i ~^ j 
which measure the movements from any municipality i 
to the municipality j and we will focus on the flows of in- 
dividuals (workers and students) commuting throughout 
the set of Sardinian municipalities by all means of trans- 
portation. This data source allows the construction of the 
Sardinian inter-municipal commuting network (SMCN) 
in which each node corresponds to a given municipality 
and the links represent the presence of a non-zero flow of 
commuters among the corresponding municipalities. 




FIG. 1: Geographical versus topologic representation of the 
SMCN: the nodes (red points) correspond to the towns, while 
the links to a flow value larger than 50 commuters between 
two towns. 

The standard mathematical representation of the re- 
sulting network is provided by the adjacency matrix A of 
elements {dij). The elements on the principal diagonal 
{an) are set equal to zero, since intra- municipal com- 
muting movements are not considered here. Off-diagonal 
terms are equal to 1 in the presence of any non-zero 
flow between i and j {i — > j or j — > i) and are equal 
to otherwise. The adjacency matrix is then symmetric 
and describes regular bi-directional displacements among 
the municipalities. The adjacency matrix contains all 
the topological information about the network but the 
dataset also provides the number of commuters attached 
to each link. It is therefore possible to go beyond the mere 



topological representation and to construct a weighted 
graph where the nodes still represent the municipal cen- 
tres but where the links are valued according to the ac- 
tual number of commuters. Analogously to the adjacency 
matrix A, we thus construct the symmetric weighted ad- 
jacency matrix W in which the elements Wij are com- 
puted as the sum of the i ^ j and j i flows between 
the corresponding municipalities (per day). The elements 
Wij are null in the case of municipalities i and j which 
do not exchange commuting traffic and by definition the 
diagonal elements are set to zero . According to the as- 
sumption of regular bi-directional movements along the 
links, the weight matrix is symmetric and the network is 
described as an undirected weighted graph. The weighted 
graph provides a richer description since it considers the 
topology along with the quantitative information on the 
dynamics occurring in the whole network. It is however 
important to stress that while the nodes correspond to 
municipalities located in the physical space, the graph 
representation does not contain any information explic- 
itly related to geographical distances and other spatial 
characteristic of the network. The definition of a net- 
work representation that correlates topological and traffic 
characteristics with the spatial properties of the SMCN 
requires more refined investigation and will be presented 
elsewhere (De Montis et al, forthcoming). 



III. THE ANALYSIS OF THE TOPOLOGICAL 
PROPERTIES OF THE SMCN 

In this section, we analyze the topology of the SMCN 
and we propose possible territorial explanations for our 
findings. This network is relatively small being charac- 
terized by = 375 vertices and E = 16, 248 edges. The 
degree k of nodes, which measures the number of links 
of each node, ranges between 8 and 279 and exhibits a 
high average value . The degree can be considered as a 
first indication of the topological centrality of a vertex: 
the importance of the corresponding municipal centre as 
attraction point for commuters in the network. A mea- 
sure of the topological distance within the network is 
given by the shortest path length defining the distance 
between nodes in terms of links to be traversed by us- 
ing the shortest possible path between the nodes. On 
the average, the shortest path length is I ~ 2.0, while 
the maximum length is just equal to 3. These values 
are small compared to the number of nodes, in agree- 
ment with a small-world behavior for which 1 typically 
scales as the logarithm of N. Further information on 
the network's topology is provided by the degree distri- 
bution P{k) defined as the probability that any given 
node has degree k. In this case it is essential to dis- 
tinguish between an 'exponential' network similar to the 
usual random graph for which P{k) decreases exponen- 
tially fast, and scale- free networks indicative of hub-like 
hierarchies for which P{k) decreases typically as a power 
law (Albert and Barabasi, 2002; Amaral et al, 2000). In 
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the present case, the degree probabihty distribution is 
skewed and is relatively peaked around a mean value of 
order 40 (Figure 2). The tail of the distribution contains 




FIG. 2: Plot of the probability distribution of connectivity. 
The red line is a lognormal fit. 

outliers with respect to a regular exponential decay, sig- 
naling the presence of municipal hubs, which exchange 
commuters with many others municipalities. However, a 
closer inspection of the distribution tail shows that the 
decay is not compatible with a heavy-tailed distribution. 
The tail has indeed a power-law behavior however with 
an exponent that ensures the convergence of the relevant 
moments of the distribution and does not imply a scale- 
free heterogeneity. In particular, the variance converges 
to a finite limit for infinitely large systems in contrast 
with scale-free networks, which have an infinite variance 
for infinite size networks. 

A further relevant topological quantity is the clustering 
coefficient defining a measure of the level of cohesiveness 
around any given node. It is expressed as the fraction of 
connected neighbors 



C{^) 



2E, 



1) 
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where E(i) is the number of links between the neighbors 
of the node i and ki(ki — l)/2 is the maximum num- 
ber of possible interconnections among the neighbors of 
the node. This quantity is defined in the interval [0, 1] 
and measures the level of local interconnectedness of the 
network. If C{i) = 0, the neighbors of the node i are 
not interconnected at all, while C{i) = 1 corresponds to 
the case where all the neighbors are interconnected. A 
large clustering thus indicates the existence of a locally 
well-connected neighborhood of nodes. It is often conve- 
nient to average C{i) over all nodes with a given degree 
k obtaining the clustering spectrum 
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where NP{k) is the total number of nodes of degree k. 
In Figure 3, on the left, we plot the clustering spectrum 



FIG. 3: On the left: Scatterplot of the clustering coefficient 
versus degree. On the right: assortativity of the SMCN show- 
ing a slight disassortative behaviour. 



of the SMCN: we observe a decaying behavior with C(fc) 
that varies from 0.8 to values lower than 0.2. This indi- 
cation of a decreasing clustering suggests that the SMCN 
behaves like other studied real networks, in that low de- 
gree municipalities belong to well interconnected working 
or studying communities while the high degree vertices, 
the urban hubs, connect otherwise disconnected regions. 
It is worth noting that this result is similar to the one 
obtained by Chowell et al (2003) for urban movement. 
There is a striking difference however since the SMCN 
is topologically very close to a random graph while the 
network studied in Chowell et al (2003) is a scale-free 
network. This similarity with a random graph is also 
confirmed by the average value of the clustering coeffi- 
cient (C) — 0.26, that results to be of the same order of 
the clustering coefficient theoretically computed for the 
case of a generalized random graph, which is 



{C)rg = 



N{kY 



(3) 



predicting the value {C)b.g ~ 0.24. 

Another important property is the degree similarity of 
the neighbors of a node. This is measured by the average 
degree of the nearest neighbors of a given node i 



kn 



(4) 



where V{i) denotes the set of neighbors of i. Analogously 
to the clustering coefficient, we can average the assorta- 
tivity over nodes with a given degree leading to 



knn{k) — 



1 



NP{k) 
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(5) 
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This spectrum measures the tendency of vertices to be 
connected with vertices with the same degree properties. 
If knn{k) increases with fc, nodes have a tendency to con- 
nect to nodes with a similar or larger degree and the 
network is said to be assortative. In contrast, if fc„„(fc) is 
decreasing with k, small degree nodes connect preferen- 
tially to hubs and the network is defined as disassortative. 
The result obtained for the SMCN is shown in Figure 3, 
on the right, and displays a fc„„ with decreasing trend. 
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This result seems to indicate that the SMCN belongs to 
the class of disassortativc mixed networks. The disas- 
sortative behavior is typical of technological and trans- 
portation network and signals the presence of a global 
hierarchy in which the hubs tend to provide the connec- 
tivity to the small degree nodes at the periphery of the 
network. 



IV. THE SMCN WEIGHTED NETWORK 
ANALYSIS 

The ODT dataset provides the number of commuters 
among municipalities, allowing the construction of the 
weighted graph representation of the SMCN. The anal- 
ysis of these commuting flows informs us on the system 
dynamical behavior and on the needs of the commuters 
in terms of regional transportation services and infras- 
tructure. This information might provide relevant indi- 
cation to researchers and decision-makers for addressing 
environmental policies and planning. 

A. Weight analysis 

The values of the weights between pairs of vertices 
ranges between 1 and Wmax = 13, 953 with an average 
value (w) « 23 much lower than Wmax ■ This fact is a sig- 
nature of a high level of heterogeneity in the SMCN, since 
a peaked distribution of weights would very unlikely give 
a maximum value that is 3 orders of magnitude larger 
than the average. Indeed, the probability distribution 
P{w) that any edge has a weight w displays a power-law 
decay P{w) ~ w~''^ with an exponent 7^ « 1.8 (Figure 
4, on the left). 
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FIG. 4: On the loft: Log-Log plot of the probability distri- 
bution of the weights (the straight line is a power law fit of 
exponent 1.8). On the right: Log-Log plot of the probability 
distribution of the strength (as a guide to the eye we also plot 
a power law fit of exponent 2). 

Another relevant quantity that characterizes the traffic 
is given by the strength of a node, defined as 

s{i)= ^ij (6) 
jev{i) 

This quantity corresponds to the total traffic commuting 
in the town i and can be considered as a measure of the 



traffic centrality of a municipality. For the SMCN, the 
strength varies between 1 and 64, 834 while its average 
value is equal to 1,010. The probabihty distribution of 
the strength P{s) that any given node has a strength s is 
fitted by a power-law regime P{s) ~ s~'^' with exponent 
7s ~ 2 (Figure 4, on the right). The strength provides an 
intuitive parameter for the description of the centrality 
of Sardinian municipalities according to the actual traffic 
handled generated and received by the municipality. 

In the first instance, the variation in the weight and 
strength values implies that commuters' flows are dis- 
tributed in a broad spectrum with a large heterogeneity. 
It is worth commenting the implications of a power-law 
behavior of the statistical distributions P{s) and P{w). 
This behavior defines the so-called heavy tailed distribu- 
tions which have a virtually unbounded variance. A pe- 
culiar fact about a distribution with a heavy tail is that 
there is a finite probability of finding vertices with weight 
or strength much larger than the average value. In other 
words, the consequence of heavy tails is that the average 
behavior of the system is not typical. The characteristic 
weight (or strength) is the one that, picking up a vertex at 
random, should be encoimtered most of the times. This 
is evident in the bell-shaped distribution observed for the 
degree of the SMCN, in which the average value is very 
close to the maximum of the distribution and represents 
the most probable value in the system. On the contrary, 
the distributions shown in Figure 4 are highly skewed. 
Vertices have an appreciable probability of having large 
strength values, yet all intermediate values are probable 
and their mean does not represent any special value for 
the distribution. The power-law behavior and the rela- 
tive exponent thus represent a quantitative measure of 
the level of heterogeneity of the network's traffic. 

Another important question concerns the distribution 
of the traffic among the different roads. All the con- 
nections could carry a similar flow or on the contrary 
one connection could dominate. A convenient measure 
of this is given by the disparity (Barthelemy et al, 2005), 
defined as 

For hubs {k 1), this quantity enables to distinguish 
situations (Figtire 5, on the right) for which all weights 
are of the same order {Y2 ~ 1/fc <C 1) from situations 
where only a small number of connections dominate {Y2 
is of order 1/n where n <^ k). 

We computed the value of Y2 for each node and we 
average these quantities for each degree in order to ob- 
tain Y2{k). The result is shown on Figure 5, on the 
left. We obtain here an average behavior of the form 
kY2ik) - k^-^ with 6 ^ 0.4. This result shows that 
the weights on the links attached to a given node are 
very heterogeneous. In other words, there are only a few 
dominant connections, while the traffic on all the other 
roads is very small. 
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Yj=0.46 



FIG. 5: On the left: Log-Log plot of the disparity versus 
connectivity for the SMCN; on the right: illustration of the 
disparity for two very different cases. When a few connections 
dominate, Y2 is of order 1; in contrast if all connections have 
the same weight, Y2 is of order 1/k. 



B. The analysis of traffic-topology correlations 

In order to inspect in more detail the relation be- 
tween the degree and the traffic, we show in Figure 6 
the strength of nodes as a function of their degree. De- 
spite the existence of inevitable fluctuations, we observe 
over a wide range of degrees a power-law behavior of the 
form s{k) ~ with an exponent (3 ~ 1.9. Indepen- 
dent weights and connectivities would give a value (3=1 
(Barrat et al, 2004) and the result obtained here reveals 
that there is strong correlation between the traffic and 
the topology. The strength of the nodes grows faster 
than their degree: the more a municipality is connected 
with other centres the much more it is able to exchange 
commuters' flows or, in other terms, the traffic per con- 
nection is not constant and increases with the number of 
connections 




FIG. 6: Average strength of the municipalities as a function 
of the degree. 

The presence of weights also alters the definition of 
other structural measures. For instance, the standard 
topological clustering coefficient used in the previous sec- 
tion needs to be augmented by the inclusion of the actual 
traffic. In order to understand the weight structure and 
their relation with topology, the weighted clustering co- 
efficient can be defined as follows (Barrat et al, 2004) 



C"'(i) 



Wih 



i (ki 



-(Xij(l{fiQ,jfi 



(8) 



where Si{ki — 1) is a normalization factor which ensures 
that C""(i) belongs to [0, 1]. As for the topological case 
it is possible to average over all nodes of same degree k 
obtaining the clustering spectrum C'^"{k). The weighted 
clustering coefficient counts for each triple formed in the 
neighborhood of the vertex i the weight of the two par- 
ticipating edges starting from i. In this way, we are not 
just considering the number of closed triangles but also 
their total relative weight with respect to the vertex' 
strength. In the case of random networks C™ = C but 
in real weighted networks, we can however face two dif- 
ferent situations. If C™ < C the topological clustering 
is generated by edges with low weight and therefore the 
cohesiveness is less important in terms of traffic proper- 
ties. On the contrary, if C™ > C we are in presence of 
a network in which the interconnected triples are more 
likely formed by edges with larger weight. We show in 
Figure 7 (on the right) an example of such a situation. 
It is clear that the interconnected triple has a major role 
in the network dynamics and organization, and that the 
clustering properties are not completely exploited by a 
simple topological analysis. The weighted clustering co- 
efficient takes into account the weighted elements by an 
appropriate mathematical form. 

In the SMCN the values of the weighted clustering co- 
efficient are larger than the corresponding topological val- 
ues > C (Figure 7, on the left) over the entire range 
of degree values. In contrast with the topological cluster- 
ing, the weighted clustering coefficient is approximately 
constant over the whole range of connectivity. This re- 




where aim is an element of the adjacency matrix and 



FIG. 7: On the left, top: weighted clustering coefficient ver- 
sus connectivity; on the left, bottom: relative difference be- 
tween the weighted and the topological clustering coefficients 
versus connectivity; on the right: large cities are well inter- 
connected by large flow links (bold edges must be considered 
with a larger weight). 

suit reveals the existence of a rich-club phenomenon in 
that important cities form a group of well interconnected 
nodes and that weight heterogeneity is enough in order 
to balance the lack of topological clustering. The rich- 
club phenomenon implies that in the SMCN the intercon- 
nected triplets are more frequently built by edges with a 
higher weight (Figure 7, on the right): in terms of local 
cohesiveness of workers' and students' commuting sys- 
tem, two different destinations available from a given city 
are more likely to be connected if the traffic flow leading 
to them is large. Similarly to the weighted clustering, we 
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need to integrate the information on weights in a sensi- 
ble definition of assortativity/disassortativity. This can 
be easily done by introducing the weighted average de- 
gree of the nearest neighbors of a given node i 



. ..... 



(9) 




Rank 


Municipal centers 


Degree k 


1 


Cagliari 


279 


2 


Nuoro 


196 


3 


Oristano 


190 


4 


Macomer 


171 


5 


Sassari 


151 


6 


Quartu Sant'Elena 


140 


7 


Assemini 


139 


8 


Selargius 


127 


9 


Villacedro 


125 


10 


Ottana 


124 



TABLE I: Ranking of Sardinian municipalities by their degree 
k. 



FIG. 8: On the left, top: weighted assortativity; on the 
left, bottom: relative difference between the weighted and the 
topological assortativity; on the right: a difference between 
the weighted and the topological assortativity indicates that 
large cities are connected by large flow links. 

In this case, we perform a local weighted average of 
the nearest neighbor degree according to the normalized 
weight of the connecting edges. This definition implies 
that k^^ > knn if the edges with the larger weights 
are pointing to the neighbors with larger degree and 
^nn < ^nn in the oppositc case. The fc^„ thus mea- 
sures effective affinity to connect with high or low degree 
neighbors according to the magnitude of the actual in- 
teractions. In the SMCN the weighted average degree 
of the nearest neighbors displays a slightly increasing 
pattern, which is a signature of assortativity (Figure 8). 
This result is in sharp contrast with the behavior of the 
topological disassortativity and indicates an affinity of 
high-degree municipalities with other large centers with 
which they exchange high number of commuters: hubs 
have more frequently large traffic flows with municipal- 
ities with the same or higher-order degree. This phe- 
nomenon is illustrated on Figure 8, on the right: large 
weights connect hubs while the connection with nodes of 
smaller degree has smaller weights. 

V. DISCUSSION: THE SMCN AND ITS 
RELATION WITH REAL ENVIRONMENTAL 
SYSTEMS 

A. Topology and environmental properties 

The ranking connected to the degree (Table 1), an in- 
dex of topological centrality, confirms the existence of ad- 
ministrative and socio-economic hierarchies: among the 
hubs, Cagliari is the most populated and regional capi- 
tal town, Nuoro, Oristano and Sassari are province cap- 
ital towns, Macomer and Quartu Sant'Elena are emerg- 
ing productive and residential centres. While the SMCN 



network is not scale free, the analysis of the degree prob- 
ability distribution still reveals the presence of a rela- 
tively large number of well-connected towns that attract 
commuters from a fairly high number of satellite cen- 
ters. The behavior of the clustering coefficient versus 
the degree indicates that for A: < 50 the clustering is 
relatively high and is a signature of the fact that small 
towns (small k) form very well connected clusters of mu- 
nicipalities exchanging students and workers. On the 
contrary, hub-towns {k > 50) have less cohesive neigh- 
borood since they have commuting flows with several dif- 
ferent local clusters of municipalities not interconnected 
among them. This kind of star-like topology with large 
municipalities acting as hub for several smaller groups 
of municipalities is confirmed by the analysis of the av- 
erage nearest neighbor degree spectrum suggesting that 
the SMCN belongs to the class of dissassortative mixed 
networks. In this case, high degree hubs link preferen- 
tially with low degree municipalities. This topological 
analysis would suggest the existence of inter-municipal 
commuting districts of small towns, which pivot around 
a few urban poles and could be interpreted as the sig- 
nature of a phenomenon often observed in urban and 
transportation systems when top functional rank towns 
attract commuters from small centers, behaving as urban 
poles. In such an environment, citizens of the satellite 
centers, being interested to regional-level services such 
as public administration, health and finance, prefer com- 
muting to the pole rather than to a peer level town. The 
heterogeneity detected in the pattern of municipal cen- 
ters could be seen as functional to an efficient behavior 
of the whole network, since each urban hub exchange 
commuters, and thus goods, services and revenues, with 
several small size towns. Turning to the geographical 
displacement of the hub centers, we can note that as ex- 
pected they are located far away from each other. The 
distance between Cagliari and Sassari is equal to 210 Km, 
while from Cagliari and Nuoro to 160 Km, with respect 
to a total length of the island equal to 230 Km. 
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B. Traffic and environment in tlie weighted 
characterization of the SMCN 

The heterogeneity of the weights, as we saw in the pre- 
vious chapters, is very high. This imphes that the hnks 
are not all equivalent and thus warns against any inter- 
pretation based on the sole topology. Indeed, the analysis 
of the weighted network introduces new insights in the 
interpretation elaborated on the result of the topological 
properties under a variety of perspectives. 



Rank 


Pairs of connected municipal centers 


Weight w 


1 


Cagliari-Sassari 


13,953 


2 


Sassari-Olbia 


7,246 


3 


Cagliari-Assemini 


4,226 


4 


Porto Torres-Sassari 


3,993 


5 


Cagliari- Capoterra 


3,731 



TABLE II: Ranking of pairs of Sardinian municipal centres 
by the weight of their connections. 



Since the weights are distributed in a very heteroge- 
neous pattern, the corresponding commuters' flows dis- 
play a wide range of values, which reflects the existence 
of a territorial hierarchy embedded and supported by the 
underlined transportation network. Precise indications 
emerge about dominant relations between some of the 
major poles of the island (Table 2): Cagliari, Sassari, 
Olbia, Assemini, Porto Torres and Capoterra. 

In particular, the analysis of the strength probability 
distribution shows that municipalities identified as topo- 
logical poles behave also as attractors of traffic flows. 
The strength is broadly distributed along the nodes, re- 
vealing a strong hierarchical dominance of the dynamic 
hubs over other smaller level centers. The analysis of the 



Rank 


Almiicipal cciitt-rs 


StiTugth s 


1 


Cagliari 


64,834 


2 


Sassari 


21,437 


3 


Quaxtu Sant'Elena 


18,431 


4 


Oristano 


12,130 


5 


Sclargius 


10,084 


6 


Assemini 


7,915 


7 


Porto Torres 


6,886 


8 


Nuoro 


6,834 


9 


Carbonia 


6,616 


10 


Iglcsias 


6,479 



TABLE III: Ranking of Sardinian municipalities by their 
strength. 



relation between strength and degree shows that traffic 

grows superlinearly with the degree confirming the fact 
that the more the topological hub-towns are connected 



and the larger their traffic centralities (ie. their capac- 
ity to attract Sardinian commuting students and work- 
ers). The ranking of Sardinian municipalities according 
to their strength is reported in Table 3. 

From the behavior obtained for the disparity (Fig. 
5), we see that only few connections emerging from the 
hub-towns carry a large amount of traffic, indicating a 
strong hierarchy and heterogeneity also in the local struc- 
ture of the flows connecting large municipalities of the 
SMCN. This mirrors the underlying road network struc- 
ture, which comprises a very small set of main highways 
(state roads), a large set of medium-sized roads (provin- 
cial roads) and a very large set of local roads (munici- 
pal roads). In particular, the link Cagliari-Sassari corre- 
sponds to state road number 131, named 'Carlo Felice' 
that is the highway crossing longitudinally the island and 
can be considered as the backbone of the SMCN struc- 
ture to which the lower hiearchical levels transportation 
structures connect. 

The use of the traffic information provides further in- 
sight on the structure and hierarchy of the network. In 
particular the difference detected between weighted and 
topological clustering coefficients and weighted versus 
topological assortativity indicates the presence of a dis- 
tribution of traffic flows that might be indicative of what 
sometimes is referred to as the rich-club phenomenon. 
Indeed, it is possible to observe that large traffic flows 
are characterizing the commuter traffic among munici- 
palities with high centrality ranks as measured in terms 
of degree and strength. The accumulation of weight on 
these centers attenuates the decaying of the clustering 
and average nearest neighbour spectrum that appears to 
invert their behaviour by using the weighted definitions. 
This is the sign of the existence of a backbone of com- 
muting patterns among large municipalities which are 
all interconnected, while the flows toward peripheral mu- 
nicipalities are indeed of small size. According to this 
picture, small municipalities are mostly connected each 
other through small traffic links, which correspond to 
the underlying second or third order roads. This com- 
pletes the topological information and offers a clear pic- 
ture of the ordering principles of the SMCN, in which 
a large number of small-size cities tend to become satel- 
lites of higher rank municipalities, shaping an overall net- 
work structure widely punctuated with star-like subsys- 
tems pivoting around important urban poles connected 
through dominant roads and the corresponding traffic 
flows. In this respect it appears that the network analysis 
presented here is able to provide an extremely detailed 
and quantitative characterization of the structure of the 
inter-urban traffic system. 



C. Relating topology and traffic to socio-economic 
phenomena 

In this section, we provide a discussion of the emergent 
topological and traffic properties in relations to the de- 
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mographic and economic properties of the towns of the 
SMCN. In the previous sections, we observe that the val- 




FIG. 9: On the left: map of the degree k; on the right: map 
of the strength s for each Sardinian municipality (Scale: the 
darker the shade, the higher the value). 

ues of degree k and of the strength s, as indexes of topo- 
logical and weighted centrality, confirm quantitatively 
historical structural hierarchies among Sardinian towns. 
The heterogeneity of the strength and the degree can be 
visualized as shown in the thematic map representations 
of Figure 9. 

We explore now in quantitative terms the relation be- 
tween the network properties and the environmental and 
economical indicators. One goal is to inspect whether 
the commuting network hubs correspond to economically 
relevant towns. As a measure of social and economic cen- 
trality of each town, we consider two variables: the total 
resident population (pop) in 1991 (Istat, 1991b) and the 
average monthly income (mmi) . The average monthly in- 
come, an index of urban aggregated wealth, is defined as 
the product of the average monthly income per worker 
times the number of workers (Carcangiu et al, 1993). 
A map representation using these indicators is shown in 
Figure 10. 

In Figure 11, we plot the variable pop as a function of 
the degree k and the strength s and we observe a clear 
positive correlation. The association between the popu- 
lation and these variables has a power-law behavior with 
two different exponents related by a factor 1.9. This is 
expected since we have already observed that traffic and 
degree are related by an empirical law s k^'^. This implies 
that if the population has a scaling pop s°, the simple sub- 
stitution of variable yields pop k^'^"". This is confirmed 
in the analysis where the slopes of the two curves are 
related by a factor close to two. It is interesting to note 
that the population scales almost linearly with the traffic 
which suggest that a constant fraction of the population 
of every city commutes to another town, irrespectively of 




FIG. 10: On the left: map of the resident population pop; on 
the right: map of the municipal monthly income mmi for each 
Sardinian municipality in 1991 (Scale: the darker the shade, 
the higher the value). 



the size of the city. 




FIG. 11: On the left: Log-Log plot of population over the 
degree k, slope coefficient equal to 1.70; on the right: Log-Log 
plot of population over the strength s, slope coefficient equal 
to 0.90. 

In Figure 12, the monthly municipal income (mmi) 
is plotted versus the degree k and the strength s. A 
positive correlation is signaled by a linear behavior on 
a log-Log scale which indicates the presence of a power 
law association. Also in this case, the exponents of the 
two relations confirm the general result of a quadratic 
relation between strength and degree. 

The results of this section confirm that the central- 
ity measures (strength and degree) display trends consis- 
tent with the behavior of social and economic indicators 
such as population and income. In this sense, the higher 
the network centrality of a node, the higher is its demo- 
graphic and economic size. The aforementioned territo- 
rial size measures, pop and mmi grow much faster with 
the degree k than with the strength s. We can therefore 
identify a group of variables related to population (s, pop, 
mmi) that are correlated by association laws with a be- 
havior very close to the linear one. Their relation with 
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FIG. 12: On the left: Log-Log plot of mmi versus the degree 
k, the slope coefficient is equal to 1.80; on the right: Log- 
Log plot of mmi versus the strength s, the slope coefficient is 
equal to 0.90. 



the degree is however highly non-linear indicating that 
these variables grows much faster than the connectivity 
of municipalities. Indeed, inverting the relation between 
degree and traffic (or population) we have k yG. The 
degree of a municipality is both an indicator of the level 
of economic and social exchange with other municipal- 
ity and the transport infrastructure. The above relation 
implies that the economic and transportation infrastruc- 
tures seem to grow at a pace slower than population in- 
dicators. It is not clear if this is a peculiar character- 
istic of the Sardinia region, due to political and urban 
planning decision, or a general property of population 
flows indicating a kind of adaptive equilibrium of a soci- 
etal system. It is interesting to note that the worldwide 
airport network exhibits association trends among the 
network infrastructure and the population size very sim- 
ilar to those observed in this study (Barrat et al, 2004), 
suggesting a possible level of universality. The detailed 
study of other commuting and transportation networks 
is however needed in order to gain further insight on this 
issue. 



VI. CONCLUSIONS 

In this paper, the inter-urban commuting system of the 
island of Sardinia has been represented as a weighted net- 
work. This representation allows the quantitative char- 
acterization of the capacity of attraction of each urban 
center on workers and students. By comparing the dif- 
ferent levels of hierarchies emerging in the network with 
the pattern of demographic and economic poles we shed 
light on several structural and organizational properties 
of the inter-municipal traffic system. 

At first sight, the SMCN is similar to a small- world 
random (Watts and Strogatz, 1998). The clustering or- 
ganization of the SMCN however diverges from the usual 
random paradigm and similarly to other technological 
networks, such as the airline network, suggests a hi- 
erarchy of nodes where small municipalities are locally 
densely interconnected, while large municipalities (re- 
gional hub) provide the long range connectivity for dif- 
ferent areas of the islands otherwise disconnected. This 



picture is reinforced by the evidence of a disassortative 
mixed network where the hub municipalities are found to 
act as star-like vertices connecting on average municipal- 
ities with low degree and low centrality ranking. 

The weighted representation of the network allowed 
by the inclusion of traffic data complements the result- 
ing picture with important information relating the com- 
muter traffic to the available topological connectivity. We 
find that both the weights and the strengths are very 
broadly distributed, conffiming the necessity of including 
those measures in a realistic description of the SMCN. In 
particular, the probability distributions of weights and 
strengths (commuter traffic handled by the municipal- 
ity) display a power-law regime over a wide range of de- 
gree values: the SMCN behaves as a scale-free weighted 
network. In addition, the traffic is strongly correlated 
to the connectivity of the municipality and moreover, as 
shown by the disparity, it is largely concentrated on a 
very few links. This indicates that the more a town is 
connected, the much more it attracts commuters along 
a small number of dominant roads. Our analysis of 
weighted clustering suggests that important towns form 
a group of well-interconnected centers. It also shows that 
in the SMCN the main fraction of traffic exchanges oc- 
curs among the hubs of the network: hub centers accu- 
mulate their commuter traffic to same or higher order 
municipalities. In terms of environmental system anal- 
ysis, a large number of small-size centers are commut- 
ing satellites of higher rank municipalities. The SMCN 
appears as widely punctuated with star-like subsystems 
pivoting around important urban poles connected prefer- 
entially through high traffic roads. In addition, the fact 
that the weighted clustering coefficient is large for the 
whole range of connectivity suggests that on the local 
scale there are micro-commuting basins, rationalized as 
inter-municipal districts of highly clustered small towns, 
which pivot around a few urban poles belonging to macro 
commuting basins. 

Interestingly enough, there are striking similarities for 
the weight distribution and the clustering obtained in 
the present study with those obtained for the networks 
of population movements at the urban level (Chowell et 
al, 2003): this evidence suggests the existence of a set of 
basic mechanisms leading to common features in trans- 
portation systems at various granularity levels. In Table 
4, the main characteristics of different transportations 
networks at different scales are shown. In particular, it 
seems that while the topological properties can vary, the 
broadness of the traffic distribution is a common feature 
to all these networks. Another common feature is the 
existence of correlation between the traffic and the topol- 
ogy. 

Finally, the existence of the municipality hierarchy 
that we have uncovered with the help of a weighted net- 
work analysis is confirmed by comparing demographic 
and economic urban polarization in Sardinia. Popula- 
tion and wealth, as proxies of the endowment of local 
resources, display a positive correlation with topological 
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Network | P{k) | P(s) | s ~ fc^ [^2 ^ fc~^ 
Global: WAN l Heavy tail I Broad I /j = l.s l 6 = 1.0 
Inter-Cities Light tail Broad [3 = 1.9 9 = 0.4 
Intra- urban Heavy tail Broad /3 = 1.0 n/a 

TABLE IV; Comparison of properties of transportation net- 
works. 



and dynamic centrality, strengthening the value of the 
network approach as a general tool for the analysis of 
inter-urban transportation systems. 

VII. OUTLOOK 

The results of this investigation open further questions 
that are currently under study. The indication of net- 
work centrality provides a powerful measure of similarity 
of the municipalities in attracting workers and students 
daily, more work however is needed to tackle the prob- 
lem of the precise identification of the communities of 
workers/students commuting in Sardinia and their ge- 
ographical basins. These studies may provide insights 
for bottom-up processes of decision-making and planning 
based on the definition of the existing emerging behaviors 
of the citizens. Finally, the SMCN is obviously not inde- 
pendent from geography and space. It refers to a phys- 
ical network, mostly based on road systems, which dis- 
play precise geographical characteristics. Further studies 
which include spatial aspects of the network are currently 
under progress. This forthcoming analysis will hopefully 
provide decision-makers and planners with relevant indi- 
cations on the global effect of new infrastructure. 
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